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SUMMARY 

Recent research activities and accomplishments at Lewis Research Center 
on composite mechanics for engine structures are reviewed in summary form. 

The activities mainly focused on developing procedures for the computational 
simulation of composite intrinsic and structural behavior. The computational 
simulation encompasses all aspects of composite mechanics, advanced three- 
dimensional finite-element methods, damage tolerance, composite structural and 
dynamic response, and structural tailoring and optimization. 

INTRODUCTION 

The application of composites in engine structures offers a multitude of 
advantages to the use of metals. These advantages include high stiffness, 
light weight, improved fatigue, higher damping, reduced number of fabricated 
parts, reduced number of joints, and minimum or no final machining. Disadvan- 
tages include a larger number of thermal and mechanical properties, low impact 
resistance, low damage tolerance, limited engineering data base, and more com- 
plex analyses. The complex analyses are conventionally referred to as compos- 
ite mechanics. Continuing research activities at Lewis Research Center have 
led to significant developments in composite mechanics for application to 
engine structures. The objective of this report is to describe recent research 
activities and accomplishments in composite mechanics for fiber-reinforced, 
polymer-matrix composites. 

These activities included research in all aspects of composite mechanics, 
advanced finite-element methods, composite fracture, and structural tailoring. 


Each research activity focused on <1) capturing the physics of composite behav- 
ior at the level of interest, (2) developing an appropriate mathematical model, 
and (3) developing the respective computer code for the computational simula- 
tion. Specifically, these items include (1) integrated composite mechanics 
analysis and the integrated composite analyzer (ICAN) computer code, (2) sim- 
plified composite mechanics for strength, (3) three-dimensional finite-element 
for micromechanics, (4) composite fracture toughness and progressive fracture, 
(5) local interlaminar damage, (6) passive damping, and (7) structural tailor- 
ing. Each of these is described in summary form with illustrative examples 
and brief discussions on implications for practical applications. Symbols 
used in the report are defined in the appendix. Relevant references are cited 
for additional details and for comparison with available data. 

INTEGRATED COMPOSITE MECHANICS 

The development of composite mechanics encompasses all aspects from 
micromechanics to laminate theory. Parts of the composite mechanics also 
include stress concentrations due to open holes, free-edge interlaminar 
stresses, laminate failure, and hygrothermal environmental effects. One way 
to integrate all these aspects of composite mechanics is to include them in a 
modular, open-ended, user-friendly computer code. A computer code of this 
type, identified as ICAN for Integrated Composite Analyzer, has been developed. 
Both stand-alone and portable, ICAN is based on constituent material proper- 
ties which are available in a resident data bank. Additional features of ICAN 
are its capability to handle hybrids, both intraply and interply, and the in- 
clusion of the interply layer as a distinct matrix layer. 

The integrated analysis capability of ICAN is depicted schematically in 
Figure 1. A flowchart and a portion of the input data are shown in Figure 2. 
This version of ICAN has been documented (theory and users manual) in 
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References 1 and 2 and is available through COSMIC (contact COSMIC, The Univer- 
sity of Georgia, Athens, GA 30602, concerning the availability of this pro- 
gram). Representative results predicted by ICAN are shown in Figure 3, and 
comparisons with measured data are given in Table I (Ref. 3). The computer 
code ICAN is available to perform the various composite mechanics that may be 
required for the analysis of composite structures for propulsion systems. 

SIMPLIFIED MICROMECHANICS FOR STRENGTH 

The development of composite mechanics has often led to derivations of 
simplified sets of equations for predicting specific composite properties. 

These sets of equations are useful because all the participating variables, 
their relationships to each other, and their significance and contribution to 
the specific composite property are readily observed. A set of this type was 
developed for predicting ply uniaxial strengths by using respective constituent 
properties. A part of this set is summarized in Figure 4. The various 
strengths are readily identifiable. The required inputs to the equations and 
their respective outputs are summarized in the block diagram of Figure 4. As 
shown, fabrication process variables (fiber and void volume ratios) and envi- 
ronmental effects (temperature and moisture) are accounted for in these equa- 
tions. This set of equations is described in detail in Reference 4, where 
numerical examples and comparisons with measured data are also included. 

A similar set of simplified equations was developed for ply microstresses. 
This set includes all the equations required to predict and assess the micro- 
stresses (stress in the matrix, fiber, and Interface) when the ply stresses 
are known or have been determined by using laminate theory. The concept, pro- 
cedure, and locations are illustrated schematically in Figure 5. The subset 
of equations for predicting the microstresses due to ply transverse stress are 
summarized in Figure 6. Results predicted by using these equations when the 
applied ply transverse stress is equal to transverse (1) tensile and 
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(2) compressive strength are summarized in Figure 7. The remaining parts of 
the set for predicting ply microstresses due to other ply stresses and due to 
temperature and moisture are described in detail in Reference 5. In summary, 
the development of composite mechanics also includes simplified equations 
which describe complex composite behavior in relatively simple form. Each set 
of these equations is a module in ICAN. 

THREE-DIMENSIONAL FINITE-ELEMENT MODELING 

Parallel developments in composite mechanics include three-dimensional 
finite-element modeling of composite behavior at the micromechanics scale. 

There are several reasons for this type of modeling, two of which are (1) vali- 
dation of simplified equations and (2) more detailed evaluation of micro- 
mechanistic behavior. Results predicted from a specific three-dimensional 
finite-element model to validate simplified micromechanics equations are sum- 
marized in Figure 8, where these results are compared with those predicted by 
the respective simplified equations. Similar comparisons for other properties 
are summarized in Reference 6. The modeling procedure is described in detail 
in Reference 7. Two of the properties plotted in Figure 8 (shear modulus 
Gfi.23 and Poisson's ratio v) are rather difficult to predict and very diffi- 
cult (if not impossible) to measure. In this respect, three-dimensional 
finite-element modeling is the only practical means available. 

The type of three-dimensional finite-element model shown in Figure 8, 
referred to as superelement, is readily adaptable to investigating the effects 
of interfacial bonds, disbonds, fiber breaks, voids, and even matrix cracks on 
the microstress distribution. In addition, stresses due to temperature and 
moisture can also be studied. It suffices to say that composite behavior at 
the micromechanics scale can be studied and described in detail by using 
three-dimensional finite-element modeling and analysis methods. 
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Another form of three-dimensional finite-element modeling is used to eval- 
uate the interlaminar stress field near free edges of angleplied laminates. 

This form is usually referred to as progressive finite-element substructur- 
ing. The concept and models are schematically illustrated in Figure 9 with 
typical results for a specific laminate. The details of the procedure are 
described in Reference 8, which includes similar results for other laminates. 

It is interesting to note that interlaminar stress fields near free edges may 
initiate fracture, especially under cyclic loads, which predominate in propul- 
sion structural systems. 

COMPOSITE FRACTURE TOUGHNESS AND PROGRESSIVE FRACTURE 

In addition to other design requirements, composite structures are 
designed to be damage tolerant. Damage tolerance usually refers to the 
ability of the structural component to sustain, without repairs, a specified 
level of preexisting inadvertent damage (defects) for the designed load condi- 
tions and service life. The ability of the material to be damage tolerant is 
usually measured by its fracture toughness, which is characterized by two pa- 
rameters: (1) a defect size (crack length) and (2) a stress-intensity field 

which will rapidly propagate the defect to fracture. The fracture toughness 
of a material is generally determined experimentally. One of the convenient 
experimental methods to determine the fracture toughness of a material is the 
strain energy release rate (SERR). 

Recent research at Lewis led to the development of a versatile computa- 
tional procedure to determine the interlaminar SERR of composite structures. 
This procedure, described in detail in References 9 and 10, is summarized in 
Figure 10. A loading schematic is shown in Figure 11, and representative 
results are shown in Figure 12, where SERR ' s for each fracture mode are plotted 
as functions of crack length for different angleplied composite laminates. 
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The computational procedure can be readily used to assess composite critical 
parameters such as fiber volume ratio, as shown in Figure 13. This type of 
computational procedure provides an indirect but representative method for 
evaluating the damage tolerance of candidate structural design concepts with 
available advanced composite materials and with emerging ones which are only 
available in small laboratory quantities. 

A directly related and more fundamental research activity is the computa- 
tional simulation of progressive fracture in composite structures. This proce- 
dure tracks the participating fracture modes within the different scales 
(micromechanics and constituents, macromechanics, ply, laminate, local, and 
global) in a composite structure. The procedure is described in 
Reference 11. Some representati ve results at select scales are (1) detailed 
fracture modes at the micromechanics level (Table II), (2) intralaminar shear 
stresses at the ply level (Fig. 14), (3) predominate fracture modes at the lo- 
cal level (Fig. 15), and (4) photomicrographs showing attendant fracture sur- 
face characteristics (Fig. 16). It suffices to say that computational proce- 
dures can be developed to simulate progressive fracture in propulsion compos- 
ite structures by tracking participating fracture modes at the various scales 
in which they occur. 

LOCAL DAMAGE EFFECTS ON GLOBAL STRUCTURAL RESPONSE 

Damage tolerance of composites (especially that are due to impact) can 
also be evaluated by comparing the global structural response of undamaged and 
damaged composite structures. The type of damage induced is described by the 
user, depending on the specific design requirements. For example, the inter- 
laminar damage tolerance of composite cantilevers is shown in Figure 17 in 
terms of two global variables: (1) displacement and (2) frequencies (Ref. 12). 

The global variable is plotted as a function of extent of the delamination, as 
shown in the schematics. The results shown in Figure 17 are interesting 
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because they demonstrate that considerable relative damage must be present 
prior to significant changes (greater than 15 percent) in global structural 
response. 

One practical implication of these results is that field inspection meth- 
ods based on global structural response will not be generally successful in 
detecting relatively small changes (about 10 percent) in delamination-type 
damage. Higher vibration modes need to be examined. More precisely, the half 
wavelength of the vibration mode should approach the damage length (mode 4 and 
greater in Fig. 17(c)). Composite mechanics computational simulation can be 
used to evaluate the effect of local damage on global structural response and 
also to assess the sensitivity of field inspection nondestructive-evaluation 
techniques . 

PASSIVE DAMPING EFFECTS 

Applications of composites to engine fan blades or advanced turboprops 
subject the composites to a broad range of excitation (forced vibration) 
sources and subsequent fatigue. The conventional practice is to design the 
composite blade or turboprop to avoid the frequencies of the excitation 
sources, at least in the operating range. Since composite blades and turbo- 
props are made by stacking and bonding multiple layers, it is natural to con- 
sider passive damping to minimize or dampen the effects of the anticipated and 
unanticipated excitation sources. The passive damping can be implemented in 
terms of either interplied or constrained adhesive layers. 

The effects of passive damping on vibration modes can be readily eval- 
uated, and thereby desirable material damping characteristics can be identi- 
fied by using composite mechanics and forced-response finite-element 
analyses. The computational procedure is described in detail in Reference 13. 

A composite turboprop is illustrated in Figure 18, and some representative 
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results with possible applications to composite turboprops are shown in 
Figure 19. As shown in Figure 19, the amount of damping of the vibration mode 
amplitude can be determined, and the adhesive layer material characteristics 
necessary to provide this damping can be identified as functions of modulus 
and shear stress. 

STRUCTURAL TAILORING OF COMPOSITE BLADES 

Composite blade structural design requires the integration of several rel 
evant disciplines, including (1) composite mechanics, (2) structural analysis, 
(3) fatigue, (4) flutter, (5> impact, (6) thermal stress analysis, and (7) 
overall costs, which include material, fabrication, operation, maintenance, 
and profits. Each of the disciplines is usually handled by a different indi- 
vidual or group specializing in that discipline. The discipline task is per- 
formed sequentially by using appropriate inputs from the other disciplines. 

The process usually requires several time-consuming iterations to obtain a 
satisfactory design. An effective alternative is to integrate the various par 
ticipating disciplines and the specified design requirements into a formal 
structural optimization computer code. 

Research conducted by Lewis Research Center during the recent past led to 
the development of the computer code STAEBL (Structural Tailoring of Engine 
Blades), as described in References 14 to 16. The essential features/modules 
of STAEBL are summarized schematically in Figure 20. Additional features and 
a schematic of the model produced by the dedicated finite-element module in 
STAEBL are shown in Figure 21. Typical results obtained are shown in 
Table III. STAEBL, proven to be very effective for the design of engine 
blades, is presently being extended for turboprops which have complex internal 
configurations (Fig. 18). In addition, STAEBL has been extended to include 
simultaneous tailoring for aerodynamic and structural considerations, and it 
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has also been used to obtain designs in several situations where the conven- 
tional approach did not succeed. In all these situations, professional and 
computer times required were substantially (in orders of magnitude) reduced. 
The multi-disciplinary design process of complex engine composite structures 
has been computationally simulated and integrated by using structural tailor- 
ing and optimization concepts and methods. 

CONCLUDING REMARKS 

Continuing research at Lewis Research Center in composite mechanics for 
engine structures has led to significant developments in all aspects of compos 
ite mechanics and has resulted in computer codes for respective computational 
simulations. These include (1) integrated composite mechanics analyzer 
(ICAN), (2) simplified composite micromechanics for strength, (3) three- 
dimensional finite-element modeling/superelement for composite micromechanics, 
(4) composite fracture toughness and progressive fracture, (5) local interlami 
nar damage effects on global structural response, (6) passive damping effects, 
and (7) structural tailoring of composite blades (STAEBL). Each of these are 
briefly reviewed and summarized with representative examples. The brief summa 
ries and the illustrative examples included demonstrate that the structural 
response of complex engine components can be computationally simulated at all 
levels of composite mechanics. In addition, the multidisciplinary design proc 
ess of composite blades can be integrated into structural tailoring/optimiza- 
tion codes which result in significant professional and computer time savings. 
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APPENDIX - SYMBOLS 


a 

d 

E 

G 

G T 

k 

M 

P 

S 

T 

t 

u.v.w 

X r 

x.y.z 

r 

V 

a 

1 , 2,3 

I, II, III 


crack length 
diameter 
elastic modulus 
shear modulus 

total strain energy release rate 

volume ratio 

moisture 

load 

strength 

temperature 

thickness 

displacements 

free-edge distance 

structural reference axes 

ply angle orientation 

Poisson's ratio 

stress 

ply material axes 

opening, shearing, and tearing fracture modes 


Subscripts: 

B beam 

C compression 

c composite property 

F flexural 

f fiber 

i 

a 


10 


interface 

ply 



m 


matrix 


n normal to interface 

S shear 

s symmetric, shear tangential to interface (Fig. 7) 

T tension 

u use 

v void 

x,y,z structural reference axes properties 

1,2,3 ply material axes properties 

Superscripts : 

A microstress in intermatrix subregion A 

B microstress in interfiber subregion B 

C microstress in the interface in intermatrix subregion C 


11 



REFERENCES 


1. Murthy, P.L.N. and Chamis, C. C., "Integrated Composite Analyser 
(ICAN). Users and Programmers Manual," NASA TP-2515, 1986. 

2. Murthy, P.L.N. and Chamis, C.C., "ICAN Integrated Composite Analyzer," 
Journal of Composites Technology and Research , Vol . 8., No. 1, Spring 
1986, pp. 8-17. 

3. Ginty, C.A. and Chamis, C.C., "Select Fiber Composites for Space Applica- 
tions: A Mechanistic Assessment," Technology Vectors , SAMPE, Azusa, CA, 

1984, pp. 979-993. (NASA TM-83631). 

4. Chamis, C.C., "Simplified Composite Micromechanics Equations for 
Strength, Fracture Toughness, and Environmental Effects," SAMPE Quarterly, 
Vol. 15, No. 4, 1984, pp. 41-55. (NASA TM-83696). 

5. Chamis, C.C., "Simplified Composite Micromechanics for Predicting Micro- 
stresses," NASA TM— 87295 , 1986. 

6. Caruso, J.J. and Chamis, C.C., "Assessment of Simplified Composite 
Micromechanics Using Three-Dimensional Finite-Element Analysis," Journal 
of Composites Technology and Research , Vol. 8, No. 3, Fall 1986, pp. 

77-83. 

7. Caruso, J.J., "Application of Finite Element Substructuring to Composite 
Micromechanics," NASA TM-83729, 1984. 

8. Murthy, P.L.N. and Chamis, C.C., "A Study of Interply Layer Effects on 
the Free Edge Stress Field of Angleplied Laminates," Computers and Struc- 
tures . Vol. 20, No. 1-3, 1985, pp. 431-441. (NASA TM-86924). 

9. Murthy, P.L.N. and Chamis, C.C., "Interlaminar Fracture Toughness: 
Three-Dimensional Finite-Element Modeling for End-Notch and Mixed-Mode 
Flexure," NASA TM-87138, 1985. 

10. Murthy, P.L.N. and Chamis, C.C., "Composite Interlaminar Fracture Tough- 
ness: Three-Dimensional Finite Element Modeling for Mixed Mode I, II and 

III Fracture," NASA TM-88872, 1986. 


12 


11. Chamis, C.C., "Computational Simulation of Progressive Fracture in Fiber 
Composites," NASA TM-87341 , 1986. 

12. Grady, J.E. and Sun, C.T., "Dynamic Delamination Crack Propagation in a 
Graphite/Epoxy Laminate," Composite Materials: Fatigue and Fracture . H. T. 
Hahn, ed., ASTM STP-907, ASTM, Philadelphia, 1986, pp. 5-31. 

13. Chamis, C.C. and Williams, G.L., "Interply Layer Degradation Effects on 
Composite Structural Response," Journal of Aircraft , Vol . 22, No. 7, July 
1985, pp. 573-580. (NASA TM-83702). 

14. Platt, C.E., Pratt, T.K., and Brown, K.W., "Structural Tailoring of En- 
gine Blades (STAEBL) , " NASA CR-167949, 1982. 

15. Brown, K.W., "Structural Tailoring of Engine Blades (STAEBL) - Theoreti- 
cal Manual," NASA CR-175112, 1986. 

16. Brown, K.W., "Structural Tailoring of Engine Blades (STAEBL) - User's Man- 
ual," NASA CR— 1751 13, 1986. 


TABLE I. - REPRESENTATIVE DATA FOR VARIOUS USE TEMPERATURES T u 
PREDICTED BY INTEGRATED COMPOSITE ANALYZER 
(ICAN) COMPARED WITH EXPERIMENTAL DATA 



Longitudinal elastic modulus, ksi 

Tu = 

-300 °F 

Tu : 

= 70 °F 

Tu : 

200 °F 

ICAN 

Measured 

ICAN 

Measured 

ICAN 

Measured 

Composite 1* 

4589 


4251 

4357 

4076 

4107 

Composite 2 b 

5587 

■ 

5395 

5964 

5457 

5981 

Composite 3 C 

4440 

■ 

4114 

4300 

3948 

4200 


?7781 E-glass cloth. 
b 7576 E-glass cloth. 

c Representati ve laminate; combination of 7781 and 7576 glass. 
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TABLE II. - FRACTURE MODES OF [±0] s GRAPHITE/EPOXY 
LAMINATES PREDICTED BY COMPOSITE DURABILITY 
STRUCTURAL ANALYSIS (CODSTRAN) 

[Longitudinal tension, LT ; transverse tension, TT ; 
intraply shear, S; and intraply del ami nati on , I.] 


Ply 

ori entation f 

e, 

deg 

Fracture mode 

Unnotched 
sol id 
specimen 

Notched 
specimen , 
wi th si i t 

Notched 
specimen , 
with hole 

0 

LT 

s a 

s a 



LT 

LT 

3 

LT 

S a 

s a 


$ b 

LT 

LT 

5 

LT 

s a 

s a 


s b 

LT 

LT 

10 

LT 

$ 

S 


S b 



15 

I 

S 

S 


S 


LT 

30 

s 

i c 

i c 



s 

s 

45 

I 

I c 

I c 


s 

s 

S 




TT 

60 

TT 

I c 

I c 



TT 

TT 



S« 


75 

TT 

TT 

TT 

90 

TT 

TT 

TT 


a Initial fracture due to intraply shear in notch- 
tip zone. 

b Some intraply shear occurring near constraints 
(grips). 

c Delami nati ons occur in notch-tip zone prior to any 
intraply damage. 

^Minimal intraply shearing during fracture. 



TABLE III. - EFFECT OF PRESSURE ANO THERMAL LOADS ON OPTIMUM 

BLADE DESIGN 



Pressure 

Pressure and 



load 


thermal loads 3 


Span, percent 

Span, percent 


0 

50 

100 

0 

50 

100 

Thickness, in. 

0.47 

0.08 

0.10 

0.48 

0.08 

0.08 

Chord, in. 

3.33 

3.66 

4.22 

3.27 

3.60 

4.15 

Thickness-to- 

.14 

.02 

.02 

.14 

.02 

.02 

chord ratio 







Weight, lb. 


9.88 



9.88 


Constrai nts 







Resonance margin 







Mode I 


0.05 



0.05 


Mode II 


1.55 



1.47 


Mode III 


1.71 



1.67 


Flutter constraint 


.510 



.520 


Root stress 


.782 



.813 



a Temperature-dependent properti es . 


TO 

GLOBAL 

STRUCTURAL 

ANALYSIS 



MICROMECHANICS \ 

THEORY 

CONSTITUENTS 




COMPONENT 
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GLOBAL 

STRUCTURAL 

ANALYSIS 



NONLINEAR 

HYGROTHERMOMECHANICAL 

RELATIONSHIPS 



dM/ PLY 

✓ COMPOSITE 

MICROMECHANICS 

THEORY 


M 


MATERIAL PROPERTIES 


P = P(0, T, M) 


FIGURE 1. - ANALYSIS CAPABILITY OF INTEGRATED COMPOSITE ANALYZER (ICAN). 





(A) ICAN FLOW CHART. 


SUMMARY OF INPUT DATA 

FOUR PLY SYMMETRIC LAMINATE. ICAN SAMPLE INPUT DATA. 


CASE CONTROL DECK 

NUMBER OF LAYERS NL = A 
NUMBER OF LOADING CONDITIONS NL C = 1 
NUMBER OF MATERIAL SYSTEMS NMS = 2 


COMSAT CSANB BIDE 
T F F 


RINDV NONUDF 
F T 


LAMINATE CONFIGURATION 


PLY 

NO 

MID 

DELTAT 

DELTAM 

THETA 

T-NESS 

PLY 

I 

I 

0.000 

0.0% 

0.0 

0.010 

PLY 

2 

2 

0 .000 

0.0% 

90 . 0 

0 .005 

PLY 

3 

2 

0 .000 

0 . 0% 

90 . 0 

0.005 

PLY 

A 

1 

0.000 

0 . 0% 

0 . 0 

o . oio 


- - - 

COMPOSITE MATERIAL SYSTEMS 

- - - 




MATCRD 

MID 

PRIMARY 

VFP 

VVP 

SECONDARY 

vsc 

VFS 

vvs 

MATCRD 

1 

AS— IMLS 

0.55 

0.02 

AS— IMIS 

0.00 

0.57 

0 .03 

MATCRD 

2 

SGLAHMHS 

0.55 

0 .01 

AS — IMHS 

0. AO 

0.57 

0 .01 


LOADING CONDITIONS 

PRESCRIBED LOADS FOR THE LOAD CONDITION 1 


INPLANE LOADS 

NX 

= 

1000.0000 

LB/IN 


NY 

= 

0 .0000 

L3/IN 


NXY 

= 

0 . 0 00 0 

L3/IN 

SENDING LOADS 

MX 

= 

0 .0000 

LB. IN/IN 


MY 

r 

0 .0000 

LB. IN/IN 


MXY 

= 

0.0000 

LB. IN/IN 

TRANSVERSE LOADS 

DMX/QX 

= 

0 . 0000 

L 3/IN 


DMY/QY 

= 

0 . 0000 

L3/IN 

TRANSVERSE PRESSURE 

PU 

= 

0 . ocoo 

L3/3Q. IN 

TRANSVERSE PRESSURE 

PL 

= 

0.0000 

L3/SQ. IN 


(B) SUMMARY OF INPUT DATA. 

FIGURE 2. - INTEGRATED COMPOSITE ANALYZER (ICAN). 



RESIDUAL STRESS, ksi LONGITUDINAL MODULUS, mpsi 











(A) BLOCK DIAGRAM. 


1. LONGITUDINAL TENSION: 

2. LONGITUDINAL COMPRESSION: 

FIBER COMPRESSION: 
DELAMINATION/SHEAR: 

MICROBUCKLING: 

3. TRANSVERSE TENSION: 

A. TRANSVERSE COMPRESSION: 

5. INTRALAMINAR SHEAR: 

6. FOR VOIDS: 


2 1 

S /11T ~ k f 3 * S 6 * * fT 

1 


S illC“ k f S fC 

S fllC ! “ 10S il2S +2 - 5S mT 

c « G m 

S «11C~ / G \ 

2 

— 


1 

2 

S £2ZT ~ D " " k f’ 11 " E m /E f22 ) l s mT 

— 

— 


1 

2 

S^2 C "[ 1 -<V k r k f>< 1 - E m%2»] S mC 

— 

— 

Sins 5 "! 1 ' <Vkf ' k f* 11 ' G m /G fl2 ) ] s mS 

2l - 


<4 

s 1 

1 

S m - |l-[ 4k v /(l-kfHr] 1/2 | S m 

°o < 
.o < 

3 o° 

D O 


VOID--" 


(B) UNIAXIAL STRENGTHS, IN-PLANE. 

FIGURE 4. - SIMPLIFIED COMPOSITE MICROMECHANICS EQUATIONS FOR STRENGTH. 
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1. INTERLAMINAR SHEAR: 


S i i3s-[l-(^-k f )(l-Gm/G fl2 )]s mS 


2. SHORT BEAM: SHEAR: 


3. FLEXURAL: 


^23S 


1 -Vk^(l - G m /G f23 l 


1 - k f (1 - G m /G f23 ) 
S /13SB s ’ 1 - 5S il3S 


s I 23SB“ L5 S i23S 


b mS 


, „ 3 k f SfT 

S /11F~ c- 

"it 

S fC 


b /22F 


^ 3 [! - <VkT- k f ) a - Efn/Ef22 )] s 


mT 


1 + !lDl 

s mC 


M. FOR VOIDS: S m ~ jl-[4k y /(l -k f )w] 1/2 j S n 


11888888 If 

t j QOOOQO | l 


OOOOOO 

OOOOOO 

OOOOOO 


OOOOOO 

OOOOOO 

OOOOOO 


o o o 

o o 

O O O 


(C) UNIAXIAL STRENGTHS, THROUGH THE THICKNESS. 
FIGURE M. - CONCLUDED. 



COMPOSITE STRESSES ^c* 




PLY MICROSTRESSES REPRESENTATIVE INTRAPLY RE6I0N 

FIGURE 5. - PLY MICROSTRESSES THROUGH COMPOSITE STRESS PROGRESSIVE DECOMPOSITION. 
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. y ll2 E m 
E fll 


a mll = ( v m ' f^ )°l22 


E fi I 1 

°fll * ( u fl2 ' y il2 jr^)°i22 


°m22 = ( E m/ E 22)°i22 


a mZ2 ' {^22^122 

° f22 “( E ^ 22 / E 22)°^22 
(A) 

°m33 * l v m " v i23) ( E m l E m)°l22 
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(D) TRANSVERSE MODULUS. 


FIGURE 8. - CONCLUDED. 


(E) POISSON'S RATIO. 
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FREE-EDGE DISTANCE, X r , in. 


(C) THREE-DIMENSIONAL AND INTERPLY STRESS FIELDS AS THE FREE EDGE IS APPROACHED (20° PLY, [ ±20 ] „ AS GRAPHITE-FIBER/ 
EPOXY LAMINATE). S 

FIGURE 9. - CONCLUDED. 
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(1) DETERMINE REQUISITE PROPERTIES AT DESIRED 
CONDITIONS USING COMPOSITE MICROMECHANICS. 

(2) RUN THREE-DIMENSIONAL FINITE-ELEMENT ANALYSIS 
ON ENF (MMF) SPECIMEN FOR AN ARBITRARY LOAD. 

(3) SCALE LOAD TO MATCH INTERLAMINARY SHEAR STRESS 
AT ELEMENT NEXT TO CRACK TIP. 

(M) WITH SCALED LOAD, EXTEND CRACK AND PLOT STRAIN 
ENERGY RELEASE RATE AS A FUNCTION OF CRACK LENGTH. 

(5) SELECT CRITICAL STRAIN ENERGY RELEASE RATE 
AND CRACK LENGTH. 



FIGURE 10. - GENERAL PROCEDURE FOR PREDICTING COMPOSITE 
INTERLAMINAR FRACTURE TOUGHNESS BY USING END-NOTCH- 
FLEXURE (ENF) OR MIXED-MODE-FLEXURE (MMF) METHOD. 


z,w 
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SPAN DISTANCE. 



0 .30 .50 .70 .90 


FIBER VOLUME RATIO 

FIGURE 13. - FIBER VOLUME RATIO EFFECT ON INTERLAMINAR 
FRACTURE TOUGHNESS PARAMETERS. END-NOTCH FLEXURE (AS 
GRAPHITE-FIBER/EPOXY). 



.05 0 -.05 .1 .05 0 -.05 -.1 .1 .05 0 -.05 -.1 .05 0 -.05 

NORMALIZED INTRALAMINAR SHEAR STRESS 

(A) AT D. (B) AT C. (C) AT B. (D) AT A. 


FIGURE 14. - MAGNITUDE PROFILES FOR PLY INTRALAMINAR SHEAR STRESS 0 12 IN UNIDIRECT- 
IONAL GRAPHITE/EPOXY (T300/934) COMPOSITE WITH CENTERED THROUGH-SLIT, AT VARIOUS 
LOCATIONS SHOWN ON SKETCH. 
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FRACTURE LOAD, LB 



FIGURE 15. - PREDOMINANT FRACTURE MODES IN NOTCHED LAMINATE WITH THROUGH-HOLE. 
PRESENCE OF INITIAL OR SECONDARY MODE IS INDICATED FOR INDIVIDUAL ANGLEPLY. 
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ORIGINAL' PAGE IS 

DF PP" ^ TT ,' T I'py 



(B) [±45] s LAMINATE. 

FIGURE 16. - PHOTOMICROGRAPHS OF FRACTURE SURFACES SHOWING MATRIX HACKLES AS 
DOMINANT MICROSTRUCTURAL CHARACTERISTIC INDICATIVE OF INTRAPLY SHEARING MODE. 
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DISPLACEMENT 









MAXIMUM DISPLACEMENT/ IN. 


(1) APPROXIMATE ANALYSES 

(a) CENTRAL PROCESSING UNIT (CPU) EFFICIENT 

(b) CAPTURE IMPORTANT STRUCTURAL INTERACTIONS 



DEFORMED STRUCTURE 


(2) REFINED ANALYSES 

(a) CPU INTENSIVE 

(b) DETAILED ANALYSIS 




(B) DAMPED. DAMPING = F(E/T), WHERE 
E IS MODULUS AND T IS SHEAR STRESS. 


FIGURE 19. - COMPUTATIONAL SIMULATION OF PASSIVE DAMP- 
ING IN COMPOSITE TURBOPROPS. 


(3) ANALYSIS RECALIBRATION 

(a) MODIFY CONSTRAINT LIMITS 

(b) MODIFY APPROXIMATE ANALYSES 



OPTIMIZATION USING APPROXIMATE 
ANALYSIS 

FIGURE 20. - TAILORING PROCEDURE USING APPROXIMATE AND REFINED BLADE 
ANALYSES. STRUCTURAL TAILORING OF COMPOSITE BLADES (STAEBL) . 
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(1) DESIGN VARIATIONS REFLECTED BY SCALING 
BLADE COORDINATES 

(2) BLADE MESH IS GENERATED INSIDE STAEBL 

(3) EQUIVALENT COMPOSITE MATERIAL PROPERTIES 
OBTAINED USING LAMINATION THEORY 

<« AUTOMATIC IN-PLANE ROTATION SUPPRESSION 



FIGURE 21. - DEDICATED FINITE-ELEMENT BLADE MESH GENERATOR. 
STRUCTURAL TAILORING OF COMPOSITE BLADES (STAEBL). 
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